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Abstract 

The Gell-Mann A matrices for Lie algebra su(3) are the natural basis for the Hilbert 
space of Hermitian operators acting on the states of a three- level system(qutrit). So the 
construction of EWs for two-qutrit states by using these matrices may be an interesting 
problem. In this paper, several two-qutrit EWs are constructed based on the Gell-Mann 
matrices by using the linear programming (LP) method exactly or approximately. The 
decomposability and non-decomposability of constructed EWs are also discussed and it 
is shown that the A-diagonal EWs presented in this paper are all decomposable but by 
adding A-non-diagonal terms, one can obtain various non-decomposable EWs. 

Keywords: Entanglement Witness, Gell-Mann Matrices, Linear Programming, 
Feasible Region. 

PACs Index: 03.65. Ud 
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Entanglement is one of the interesting features of quantum systems. It is used as a physical 
resource in realization of many quantum information and quantum computation processes such 
as quantum cryptography, teleportation, dense coding and quantum key distribution [H [21 [3] 
which cannot be realized in classical physics. So a strong motivation has been raised for the 
study of entanglement detection in an operational way. Though the fields of quantum informa- 
tion and quantum computation have been mainly built on the concept of qubits, however the 
exploration of qutrits has been attracted much attention in the recent years |H|5l[6]. A qutrit, 
the simplest generalization of a qubit, is a system whose operators act on a three-dimensional 
Hilbert space. Qutrits have some interesting properties that account for their usefulness: they 
have improved the efficiency and security of many quantum information protocols, and it 
is expected that systems of entangled qutrits largely violate the non-locality via Bell tests 
[7|. Therefore, the detection of entanglement, that is, distinguishing between separable and 
entangled states, for two-qutrit systems is an important problem. 

Among the different criteria to analyze the separability of quantum states, the entanglement 
witnesses (EWs) are of special interest since it has been proved that for any entangled state 
there exists at least one EW which detects its entanglement [H [9]. The EWs are Hermitian 
operators which have non-negative expectation values over all separable states and detect some 
entangled states. Therefore, for construction of EWs, one needs to determine the minimum 
of this expectation value and to demand its non-negativity which lies in realm of the convex 
optimization problem. The linear programming (LP), as a special case of convex optimization, 
is one of the very useful approaches for constructing EWs [101 ttH [El [131 [HI [13 [IS] since it 
can be solved by using very efficient algorithms such as the simplex and interior-point methods 
( see e.g. [HI [18]). 

Inside the several problems concerning the EWs, the problem of how to construct EWs by a 
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given set of operators has a great importance. In this paper, we use the Gell-Mann A matrices 
, a basis for the Lie algebra su(3) [19], to construct EWs for two-qutrit states. Two A-diagonal 
and A-non-diagonal cases are considered and it is shown that the problem of constructing such 
EWs can be reduced to a LP problem. We show that, as far as related to the LP method, 
the presented A-diagonal EWs are all decomposable while a large number of A-non-diagonal 
EWs are non-decomposable. To show the non-decomposability of these EWs, we find a class of 
bound entangled states or used known bound entangled states that can be detected by them. 
Also, it is noted that in determining the decomposability or non-decomposability of a given 
EW, the pure product states which have zero expectation values with it play a crucial role. 

Determination of the feasible region is the main goal of any separability problem. So 
in the method of the present paper, it is tried to determine the feasible region exactly or 
approximately. When the feasible region was determined, any hyper-plane tangent to it would 
correspond to an optimal EW. In this way, we get all possible optimal EWs, including the 
pre-existent ones obtained by other methods. 

The paper is organized as follows: In Section 2, we review the basic notions and definitions 
of EWs relevant to our study and describe a general approach of constructing EWs by the 
LP method. In Section 3, we consider the construction of A-diagonal EWs and Section 4 is 
devoted to the construction of various A-non-diagonal EWs. The paper is ended with a brief 
conclusion and three appendices. 

2 EWs and LP method 
2.1 Entanglement witnesses 

First let us recall the definition of entanglement and separability. A bipartite quantum mixed 
state p e B{T-t) (the Hilbert space of bounded operators acting on H = Hd^ (S> Hrfj) is called 
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separable if it can be written as a convex combination of pure product states, that is 

p = E^'^i"J'^)<«J'^i®i«!'^)<"J'^i (2.1) 

i 

where \a[''^) are arbitrary but normahzed vectors lying in the Ti-dj, and Pi > with Xli^** ~ ^■ 
When this is not the case, p is called entangled. 

An entanglement witness is a Hermitian operator such that Tr{Wps) > for all separa- 
ble states ps and there exists at least one entangled state pe which can be detected by W, that 
is Tr(Wpe) < 0. The existence of an EW for any entangled state is a direct consequence of 
Hahn-Banach theorem [20j and the fact that the space of separable density operators is convex 
and closed. 

Based on the notion of partial transpose, the EWs are classified into two classes: decom- 
posable (d-EW) and non-decomposable (nd-EW). An EW W is called decomposable if there 
exist positive operators V, Q such that 

W = V+Q^^ (2.2) 

where Ti denotes the partial transpose with respect to partite 1 and it is non-decomposable 
if it can not be written in this form [10]. Clearly a d-EW can not detect bound entangled 
states (entangled states with positive partial transpose (PPT) ) whereas there are some bound 
entangled states which can be detected by an nd-EW. 

Usually one is interested in finding EWs W which detect entangled states in an optimal 
way in the sense that when we subtract any positive operator from W then it does not remain 
an EW anymore [21]. In other words, if there exist e > and a positive operator V such that 
W = W — eV is again an EW then we conclude that W is not optimal, otherwise it is. 

2.2 Linear programming method 

This subsection is devoted to describing linear programming (LP) and a general approach for 
constructing EWs by the LP method exactly or approximately [TTJ. 
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Let us consider a non-positive Hermitian operator of the form 

W^aJ + Y,a,Q, (2.3) 

i 

where Q. are Hermitian operators and a.'s are real parameters with > 0. In this work, 
the Q. will be considered as tesor products of Gell-Mann matrices for Lie algebra su(3). We 
will attempt to choose the real parameters such that W becomes an EW. To this aim, we 
introduce the maps 

F = Tr(Q,p,) (2.4) 

for any separable state Ps- The maps P^ map the convex set of separable states into a bounded 
convex region which will be named feasible region. The first property of an EW is that its 
expectation value over any separable state is non-negative, i.e., the condition 

Fw-.^ Tr{Wps)^a, + Y,a,R>0 

i 

is satisfied for any point (Pj, P2, ■■■) of the feasible region. In order to satisfy this condition, it 
is sufficient that the minimum value of be non-negative. We note that the quantity F^ 
achieves its minimum value for pure product states, since every separable mixed state ps can 
be written as a convex combination of pure product states, say Ps — X]jPi|7i)(7i| with Pi > 
and ^jPi = 1, whence we have 

TriWp,) = J2piTr{W\^i){ji\) > C^,„, (2.5) 

i 

with Cmin^ min Tr{W\-f){-f\) 

where Dprod denotes the set of pure product states. Here in this work we are interested in 
the EWs that their feasible regions are of simplex (or at most convex polygon) types. The 
manipulation of these EWs amounts to 

minimize Fw = + Ei a^P^ 
subject to Eiic^.P, -d.)>Q j = 1, 2, ... 



Two-qutrit EWs 7 

where c-. and d. are parameters of hyper-planes surrounding the feasible regions. So the 
problem reduces to a LP problem. On the basis of the LP method, minimum of an objective 
function Fw always occurs at the vertices of bounded feasible region. Therefore the vertices 
of the feasible region come from pure product states. 

It is necessary to distinguish between two cases: (a) exactly soluble, and (b) approximately 
soluble EWs. In the case a the boundaries (constraints on P.) come from the finite set of 
vertices arising from pure product states and construct a convex polygon, while in the case 
b the boundaries may not be hyper-planes. In this case we approximate the boundaries with 
hyper-planes and clearly some vertices of the approximated feasible region do not arise from 
pure product states. Both cases can be solved by the well-known simplex method. The simplex 
algorithm is a common algorithm used to solve an optimization problem with a polytope 
feasible region, such as a linear programming problem. Here, considering the scope of this 
paper, a complete treatment of the simplex algorithm is unnecessary; for a more complete 
treatment please refer to any LP text such as [171 [18] . 

3 A-diagonal EWs 
3.1 Exactly soluble EWs 

In this subsection we consider A-diagonal EWs which can be solved exactly by the LP method. 
Let us consider a situation that the Hermitian operator is as follows 

8 

W = a,I,^ I, + J2a.K^K (3.7) 

k=l 

where Aj.'s are the Gell-Mann matrices for Lie algebra su(3). Evidently when all of its nine 
eigenvalues are positive the above operator is positive; otherwise it may be an EW. 

Now we attempt to choose the real parameters Oj such that W becomes an EW. To this aim, 
we need to minimize the expectation value of W over all separable states and to demand its 
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non-negativity. To reduce the problem to a LP problem and to determine the feasible region, 

wc require to know the vertices, namely the extreme points, of the feasible region. Vertices of 
the feasible region come from pure product states and are listed in the table 1 where jA^^,; m) is 
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Table 1: The product vectors and coordinates of vertices for W. 



the eigenvector of A^ with eigenvalue m and P^. = Tr(Aj.(g)A^|7)(7|) in which I7) is an arbitrary 
pure product state. 

The hyper-planes passing through the two vertices (0, 0, 1, 0, 0, 0, 0, |) and (0, 0, 0, 0, 0, 0, 0, |) 
together with other vertices except the vertex (0, 0, 0, 0, 0, 0, 0, — |) have the form 

3 * 

-^(-ly.P.-l^O ; i3 = is = and e {0,1} for j^3,8. (3.8) 

To locate the position of the feasible region with respect to these hyper-planes, we note that 
the origin (P^, Pg) = (0, 0) come from the pure product state 



1 

V3 



1 
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and hence lies within the feasible region. On the other hand, for the origin we have 



iE(-ir^^.-i = -i<o- 



Hence, the feasible region lies in one side of the hyper-planes (13.81) . To make sure, it is proved 
that there exists no pure product state on the other side of these hyper-planes (for a proof, 
see appendix II). So the hyper-planes (13. 8p are boundaries of the feasible region. 

Each hyper-plane (13.81) corresponds to a positive operator or an EW. To see this, note that 
each hyper-plane (13.81) can be rewritten as 



with 1 7) a pure product state that corresponds to a vertex of the feasible region. Hence, the 
operator 



has positive expectation value over any pure product state and we conclude that it is a positive 
operator or an EW. To simplify the analysis, we consider the following cases with this in mind 
that in all cases = ig = 0. 



This case contains eight operators which are all positive. 

b- i.y^i^ , hy^i, , h 7^ h- 

The eight operators of this case are all d-EWs since they have some negative eigenvalues and 
their partial transposes with respect to the first partite which is equivalent to A. — > —A- for 
i = 2,5,7, yield one of the eight positive operators of case a. 

c- i-^ = , ^ , ig ^ i,^ , i-^ ^ , = , i^^ ^ , i-^ ^ , ^ , = i,^. 

This case contains 24 operators which are all d-EWs. For instance, consider the following one 




(3.9) 



j=i 




(3.10) 




Two-qutrit EWs 10 
It can be written as 

in which the positive operators V and Q are 

V = 3\^P,){^P,\ , Q = 3(|^,)(^J + |^3)(^3|). 

Here \ip-)^s are states that to be defined in relation 01-60 of appendix III. 

d- i-^ ^ , = ir^ , if^ = , i-^ = ^ ^ i^^ , = , i-^^ = ^ = ^ ^ i^- 

The 24 operators of this case are all d-EWs since their partial transposes are equal to one of 

the 24 d-EWs of case c. 

In summary, it was shown that among the 64 boundary hyper-planes (13. 8 p of the feasible 
region, eight of them correspond to positive operators and the others correspond to d-EWs. 
These 64 boundaries are hyper-planes by itself while as it is shown in the sequel, the remaining 
ones are not hyper-planes by itself and we have to approximate them by hyper-planes in order 
to apply the LP method. 

The EW introduced in Eq. (74) of [6] for the isotropic qutrit states and also the EWs in 
Eqs. (118) and (134) of [22j are among the EWs of case b. 

3.2 Approximately soluble EWs 

In this subsection we consider other boundaries of the feasible region for EWs of the form 
(13. 7p . We show that these are not hyper-planes by itself and we have to approximate them 
with hyper-planes. 

At first, we choose two vertices (0, 0, 1, 0, 0, 0, 0, |) and (0, 0, 0, 0, 0, 0, 0, — |) together with 
other vertices except the vertex (0, 0, 0, 0, 0, 0, 0, |). The hyper-planes passing through these 
vertices have the form 

^ [{-ir^p, + (-i)^^n + p,-p,] + \ [{-ly^p, + (-1)^.^'^ + {-ly^p, + {-ly^p,] -1 = 

(3.12) 



Two-qutrit EWs 11 

where i-^^i^^i^, G {0, 1}. In appendix II, we find pure product states that maximizes the 
left hand side of fl3.12p and show that this maximum is not zero but it is |. This implies that 
the hyper-planes (13.121) lie inside the feasible region and so are not really the boundaries of the 
feasible region. At first approximation, we approximate the corresponding real boundaries by 
hyper-planes that are parallel to the hyper-planes (13.121) and are tangent to real boundaries 
of the feasible region at pure product states giving the maximum value |. The consequence 
of the above approximation is that some vertices of the approximated feasible region do not 
come from pure product states. The EWs corresponding to these approximated boundary 
hyper-planes cirG clS follows 

W, = fl, ® 4 - |[(-1)H A, ® A, + i-ipX, ® A, + A3 ® A3 - \ ® A3] 
® \ + (-1)^5 A, ®\ + (-1)%\ ® A, + (-I)VA, ® A,]. 
Among them we study the following one in detail 

W[ = y/3®/3-^[A,®A, + A,®A, + A3®A3-A,®AJ-^[A,®A, + A,®\ + \®\ + \®A,]. 

(3.14) 

To show the decomposability of W[, we assume that it is decomposable and write 

W[ = V+Q^K (3.15) 

Also, we note that Wl vanishes on the pure product states I7) = \a) \(3) in relation (1I-4P of 
appendix II, i.e., 

Tr{Wi\j){j\) = {i\Wi\j) = (7IPI7) + (tIQ^Ht) 

= (7l^l7) + («1(/3|Q|«*)l/9) = 0, 
where \a*) is the complex conjugate of \a). This requires that 

{-f\V\-f) = and {a*mQ\a*m = 0. 

So the positive operators V and Q have to be orthogonal to all |7)'s and ® |/3)'s resppec- 
tively. As it is shown in appendix II, such V and Q actually exist and hence Wl is a d-EW. 
Similarly, we can show that EWs (13.131) are all d-EWs. 



(3.16) 



(3.17) 
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In the same way, it can be shown that the remaining boundaries are also not hyper-planes 
by itself and the EWs corresponding to hyper-planes approximating them, i.e., 

W2 = 2J3 ® J3 - I [(-1)^1 A, ® A, + (-1)^2 A, ® A, - A3 ® A3 + (-1)^4 A, ® A, 
+ (-l)*5A, ® A, + (-1)^6 A, ® A, + i-iyrX^ ® A, - A3 ® \], 

and 

Ws = fl, ® 4 - I [(-1)^1 A, ®\ + i-iy^X, ® A, - A3 ® A3 + \ ® A3] 
-|[(-1)^4A, A, + (-1)^5 A, ® A3 + (-1)^6 A, ® A, + (-l)VA, ® A,], 
are all d-EWs. 

In summary, the EWs of the form (13.71) which correspond to boundaries or approximated 
boundaries hyper-planes of the feasible region, are all d-EWs. The other EWs of the form (13. 7p 
can be written as a convex combination of a positive operator and one of the EWs coming 
from boundary hyper-planes [15]. So if all the boundaries would be hyper-planes by itself, we 
could say that A-diagonal EWs are all decomposable. Therefore, it is tempting to add some 
off-diagonal terms A. ® X. with i 7^ j to the operator of (13.71) in order to get nd-EWs which is 
the topic of the next section. 



4 A-non-diagonal EWs 

4.1 First category 

For the first ctegory of A-non-diagonal EWs, let us consider a Hermitian operator of the form 
W = aoJ3® Jg + ^J^^a.A, ® A, +a,_2A, ® A^ + a^^^A^® A, 

(4.18) 

+0^ 5A4 ® + Og^Ag ® A4 + Og^Ag (g) A7 + a^^fjA^ ® Ag 
Vertices of the feasible region coming from pure product states are listed in the table 2. 
The hyper-plane passing through the two vertices 

(0,0, 0,0, 0,0,0, ^,,0,0, 0,0, 0,0) , (0,0,l,0,0,0,0,i,0,0,0,0,0,0) 



Two-qutrit EWs 



13 



Pi P2 P3 Pa P5 Pq Pj Ps Pi,2 -^2,1 -Ri.s -^5,4 Pej Pjfi 



Product state 



±1 

±1 

±1 





±10 

±1 

±1 

±1 







,^ 

























12 



12 



12 



12 



12 



12 



12 



12 







±1 




















±1 
























±1 




















±1 



















±1 























±1 






|Ai;±l) 

|A3;±i) 
|A,;±1) 
|A,;±1) 
|A4;±i) 
|A5;±i) 
l\;±l) 
|A,;±1) 

l\;±i) 
|As;±i) 
|Ae;±l) 
|A,;±1) 



A,;±l) 



A2;±l> 



A3;±i) 



A,;±l) 



A,;±l) 



A4;±i) 
A5;±i) 
|Aa;±) 
|A,;±) 
|A3;±) 
|A4;±) 
l\;±) 
l\;±) 



Table 2: The product vectors and coordinates of vertices for W. 



along with all other vertices having a coordinate ±1, is 



-(P,±P,±P,,,±F,,,±F3±P,±n±P,,,±F3,,±P,±n±P,,,±F,,,±Fj-l = 0. (4.19) 
As the proof (see appendix II) shows, there exist pure product states that maximizes the left 



hand side of fl4.19p and the maximum value is |, greater than zero. So this hyper-plane can 
not be a real boundary of the feasible region and we need to approximate the corresponding 
real boundary with the following hyper-plane parallel to it 



^(P,±F,±F,,±F,,,±P3±F,±F,±P,,,±P,,,±F,±P,±P,,,±P,,,±P3)-^ = 0. (4.20) 
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The corresponding EW of hyper-plane (14.201) is as follows 

Wi = I/3 ® /3 - KeLi a, ® a, + a, ® a, + a, ® A, 



(4.21) 



+A4 (g)X, + A3 ® A4 + Ag O A, + A, ® AJ 
Wi is a nd-EW since it is able to detect a bound entangled state. To show this, we note that 
the density matrix 

p = i/3 ® /3 + 6(^(A, ® A, + A, ® A, + A, ® A, + A, ® A, + A, ® A, + A3 ® A3 + A, ® Ag 
+A3 ® A4 + A, ® Ag + A, ® A, + Ag ® A, + A, ® AJ + ® A3 + A^ ® AJ 

(4.22) 

where a and c are non-negative real parameters such that < c < is a PPT state. On the 
other hand, 

Tr(>Vip) = - 

^ ^'^ 4 a + 2c 

It is seen that for c > | the PPT state p can be detected by Wi. Hence, for | < c < p is 
a bound entangled state and Wi is a nd-EW. 

This category contains 2^^ EWs corresponding to 2^^ exact or approximated boundary 
hyper-planes of the feasible region. This set of EWs can be partitioned into 64 equivalence 
classes of size 2^ by the conjugation action of the exchange operator 



1 1 ^ 

n = TT/3®^3 + ^5^A,®A, 



(4.23) 



k=l 



which exchanges the two parties, and the operators T® /g, and M.^I^, z = 1, 2, 3 on it, where 
T is the transpose operator and 
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i 



(4.24) 



i 
0-1 

Any two EWs are said to be equivalent if they can be transformed into each other by the above 
five operators or their possible products. The operators M., z = 1,2,3 do not change Gell- 
Mann matrices A3 and Ag while induce the following transformations on the other Gell-Mann 
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matrices 

: Aj ^ — Aj, Aj — > A^, A4 — * Aj, A5 — > — A4, Ag — Ag, A^ — > — A^, 

M, : A, ^A„ A, ^-A„ A, -^-A„ A^ ^ -A^, \ A„ A, ^ -\, (4.25) 

M3 : A^ — » — Aj, Aj — s> — A2, A4 — » Ag, A5 ^ — A4, Ag — A7, A7 — > Ag. 

We have discussed one of the equivalence classes. As the discussion of the other classes is 
similar, we do not discuss them here. 



4.2 Second category 

For the second category of A-non-diagonal EWs, let us consider the following Hermitian oper- 
ator 

W = a„/3 ® I, + ELi ®K+ «3,8 V3\ ® A, + a,,3 ® A3. (4.26) 

Vertices of the feasible region coming from pure product states are listed in the table 3 where 
^3,8 = Tr(v^A3 ® Ag|7)(7|) and Ps.s is defined similarly. To show the efficiency of the present 
method, we try to construct nd-EWs which are capable to detect the well-known one parameter 
family of bound entangled states 

P. = |4 ® 4 + it(^i ® a, - a, ® a, + a, ® a, - a, ® a, + Ag ® Ag - a, ® A,) 
-8i(^3 ® A3 + \ ® \) - §{5 - 2b){X, ® A, - A3 ® A3), < 6 < 5, 

introduced in [23]. The Horodecki states are PPT for 1 < 6 < 4 and as it was shown 
in [251 [23j, the states are separable for 2 < 6 < 3 and bound entangled for 1 < 6 < 2 and 
3 < 6 < 4. In [2^, Doherty and et al. introduced an nd-EW which detects the Horodecki 
bound entangled states for 3 < 6 < 4. 

Let us proceed to construct an nd-EW which detects p^ for 3 < 6 < 4. To this aim, we 
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Table 3: The product vectors and coordinates of vertices for W". 

begin with the following ten vertices of the feasible region 

(1,0,0,0,0,0,0, i, 0,0), (0,-1,0,0,0,0,0, i, 0,0), (0,0,-1,0,0,0,0, i, 1,-1), 

(0,0,i,l,0,0,0,^,-i,-i), (0,0,^,0,-1,0,0, ^,-i,-i), (0,0,^,0,0, 1,0,^, I, i), 

(0, 0, 1, 0, 0, 0, -1, ^, i, i), (0, 0, 0, 0, 0, 0, 0, -| , 2, 0), (0, 0, 0, 0, 0, 0, 0, -| , 0, -2), 

(0,0,0,0,0,0,0,1,0,0). 
The hyper-plane passing through these vertices has the form 

l{Pi-P.+P,-P,+Pe-Pr) + l{Ps+Ps) + ^{Ps,s-Ps,s)-^ = 0- 

Some calculations show that the maximum value of the left hand side is |. The pure product state 
giving this maximum value lies on the boundary of the feasible region and we replace its corresponding 
point with one of the above ten vertices. Then we pass a hyper-plane through the second set of ten 
points and find pure product states maximizing it. Again, we replace the corresponding points of 
these states with suitable points of the second set. By continuing this process, we end up with ten 
vertices of the table 4. An analogous argument as above give the EW 
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Table 4: The product vectors and coordinates of vertices for W\. Here uj = exp(^) and i 



= is^3 ® ^ - 6^[2553(A, A, - A, » A,) + 5227(A, A, - A^ AJ] - i§(A, A, - A, A,) 
+fi(\ ® A3 + Ag ^ AJ - i^(A3 A, - Ag A3). 

(4.28) 

which corresponds to the hyper-plane tangent to the feasible region and parallel to the hyper-plane 
passing through the vertices of the table 4. The nd-EW W\ detects the Horodecki states for 
3.146) <b<4. This nd-EW is very similar to the nd-EW introduced in pi] . 
To construct another nd-EW which detects p^^ for 1 < 6 < 2, we replace the first three points and 
the 8th point of the table 4 with the three points 

(0, 0, -1, 0, 0, 0, 0, ^, -1, 1), (0, 0, 0, 0, 0, 0, 0, , -2, 0), (0, 0, 0, 0, 0, 0, 0, , 0, 2), 
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of the table 3 and the point (0, 0, ^,0, 0, |, 0, — —if) jq) respectively. The last point corresponds 
to a pure product state coming from interchanging the states of two parties in the pure product state 
of the 8th point of the table 4. The hyper-plane tangent to the feasible region and parallel to the 
hyper-plane passing through the new ten points, corresponds to the EW 

^'2 = ® ^ - 6^[2553(A, \ - A, A,) + 5227(A, A, - A, AJ] - i§(A, ^5 A, - A, A,) 

+m(^3 55 A3 + A3 A3) + ^{X, A3 - A3 ® A3). 

(4.29) 

The nd-EW detects the Horodecki states for 1 < 6 < f|(~ 1.854). 

It is noted that the EW (|4.29|) can be obtained from EW ()4.28p by the exchange operator 11 
defined in i.e., 

w'2 =n w'^n^ 

In fact, this category contains 13 x 2^ EWs corresponding to 5 x 2^ exact or approximated boundary 
hyper-planes of the feasible region. This set of EWs can be partitioned into 36 equivalence classes, 
20 classes of size 8 and 16 classes of size 16, by the conjugation action of exchange operator 11 and 

2 2 

the operators T /g, /g and on it, where T is the transpose operator and and 
are operators defined in (j4.24p . Any two EWs are said to be equivalent if they can be transformed 
into each other by the above four operators or their possible products. So we have studied one of the 
equivalence classes of size 16 which contains the two EWs W\ and W.^- The other classes can be 
studied similarly. 

5 Conclusion 

We have considered the EWs that can be constructed from Gell-Mann matrices for Lie algebra su(3) 
by using the Linear programming method. The general form of A-diagonal and two cases of A-non- 
diagonal EWs have been discussed in detail. It has been shown that in all considered cases the feasible 
region is a polygon approximately. In A-diagonal case, the boundaries of the polygon correspond to 
d-EWs or positive operators which in turn implies that up to our approximation, all A-diagonal EWs 
are decomposable. However, in A-non-diagonal case, a large number of boundaries of the polygon 
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correspond to nd-EWs. The decomposability of presented A-diagonal EWs has been proved by writing 
them expHcitly in decomposable form or by using the pure product states having zero expectation 
values with them. To show the non-decomposability of some EWs in A-non-diagonal case, we have 
constructed a class of bound entangled states or used known bound entangled states that can be 
detected by them. Although the method may be tedious but as it was shown in the case of Horodecki 
bound entangled states, it clarify that in principle one can construct suitable A-non-diagonal EWs 
for any two-qutrit bound entangled state. 

Since the boundaries of the feasible region in A-diagonal EWs are not all hyper-planes by itself, 
it remains as an open problem to show that A-diagonal EWs are all decomposable. 
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Appendix I: 

Lie algebra su(3) 

The real Lie algebra su(3) is the algebra of the anti-Hermitian 3x3 matrices with vanishing trace. 
The standard Hermitian basis of su(3) consists of the Gell-Mann A matrices: 
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which fulfill the condition 



and the squares of expectation values of A.'s over an arbitrary vector |a) sum to |; i.e., 

E((«|AJa))^ 



i=l 



4 
3' 



(I-l) 



Here, for completeness, we give a proof for the above equality. Let \ol) = cKoIO) + ^^i |1) + 0^2 12) be an 
arbitrary state. Then we can write 



(alAJa) = 2 Re(a*Q;j) , 

(a|A3|Q;) = |q!qP — lOiP , 

(alAJa) = 2 Im(a*a2) > 

(alA^la) = 2 Im(Q;*Q;2) , 



(alAJa) = 2 Im(a*aJ 

(a|A4|a) = 2 Re(Q*a2) 
(a|Ag|a) = 2 Re(Q*a2) 



HAsl«) = ^-^KP 



(1-2) 



whence 



8 1 
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= (|aj2 + \a,\y + 4\a,\Wa,\^ + jaj^) + i(i _ 3\a,\Y 
= {\af + |aj2)(|aj2 + \a,\^ + 4\a,\^) + ^(1 - 3\a,\^f 
noting the normalization condition 

l„ |2 1 i„ |2 1 i„ |2 -I 

\cXq\ + \ct^\ + IOjI = i 

the above statement can be written as 

(1 - |q,P)(1 + 3|a,p) + i(l - 3|a,p)2 = ^. 

Appendix II: 

Proving the equalities: 

In this appendix we give the proof of equalities for some boundary hyper-planes. 
The proof of ([3781) : 

Since the proofs are similar, we give a proof only for the following hyper-plane 

Let 1 7) = I a) 1/3) be an arbitrary pure product state. Then, by definition, we have 

8 8 

iE^.=iEHA,|a)(/3|AJ/3). 
i=i j=i 

Using the Cauchy- Schwartz inequality and P-ip . we get the result 
The proof of ([3^^ : 

Since the proofs are similar, we give the proof only for the following one 

By definition of P. 's and by noting to ()I-2p , the Ihs of this equality can be written as 

Ihs = ^[4Re(a*ajRe(/5*A) + 4Im(a*aJIm(/3*/3J + {\a,f - \a,f){m^ - 
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-^(1 - S\a,f)il - 3\P,f)] + ^[4Re(a*a,)Re(/3*/?,)4Im(a*a,)Im(/3*/3,) 

= -1 + 6|aJ|aJ|/3J|/?J cos(^„, - O'J + ^(|aj2 - - 
-^(1 - 3|aj')(l - + 3|aJ|aJ|/?J|/?J cos(^„, - + 3|aJ|aJ|A||/?,| cos(e,, - 

< -1 + 6KIKII/3JI/3J + - KP)(|/3oP - IAP) 

-^(1 - 3K|2)(i - + 3K|K||/3J|/3,| + 3KIK||AI|/3,| 



= -1 + 3|a,||/3,|(|aJ|/3J + |aJ|/3J + |a,||/3,|) + ^(|a„||/3ol + - ^(|«oll/3il " l«ill/3ol)' 



By using the Cauchy-Schwartz inequality, we have 



< -1 + 3K||/3,| + ^[(K||/?j + - (KIIAI - KII/^oD' + Kl' + - sKI'l/^.P - ^] 

By the following parametrization 



lap I = sin ^ cos ip \ = sin 9' cos (p' 
|aj = sin sin (p \Pi \ = sin ^' sin tp' 
|q;2|=cos6' 1/32 1 = cos ^' 

where 9, 9', (p, (p' G [0, f ], the above statement can be written as 



(1-3) 



3 cos 9 cos 9' + - [sin^ 9 sin^ i9' cos 2{p - p>') + cos^ 6* + cos^ 9' - b cos^ (9 cos^ 9' --]-! 

2 3 

Since the Cauchy-Schwartz inequality becomes an equality when the two vectors are parallel, we must 
have 9 = 9' and ip = ip' . Then 

lhs< 3cos2 6'-6cos''6l. 

The Ihs attains the maximum value | when cos^ ^ = i ■ 

As implied from the above consideration, the pure product states which saturate the maximum 
value I are of the form 



It) = |a) 



2 



^ cos ip ^ 

6**^1 sin 99 



2 



^ cos ip ^ 

e^^i sin(/9 



(1-4) 
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in which 6^ and 82 are arbitrary real phases. 
The proof of ([4720]) : 

We have 

S 7 

= 3[Re«ajRe(/3*A) + Im«aJIm(/3*/3J + Re«aJIm(/3*A) + Im(a> jRe(/3* A) 
+Re(a;ajRe(/3*/3J + Im(a*a,)Im(/3*/3,) + Re«a,)Im(/3*/3J + Im(a*a2)Re(/3*/3J 
+Re(a;ajRe(/3;/3J + Im(a*aJIm(/3;/3J + Re(a;aJIm(/3*/3J + Im(a;'ajRe(/3;/32^ 

+^(l«ol' - l«J')(l/?oP - lAl') + ^(1 - 3|a.n(l - 3|/3,P) - J 
= 3|aJ|aJ|/3J|/3J(cosK, -0j + sin(0,, +0;j) + 3|aJ|a,||/3J|/3J(cos(e„, -e^J 
+3|Qj|aJ|/3J|/3,|(cos(0„, + sin(0„, + +3|aJ|a,||/?J|A|(cos(e,, +sin(0,, + 0^^^ 

+^(i«oi' - i«j')(i/?oi' - m') + ^(1 - 3|a.i')(i - - 1 

15 3 

< elaJlaJdaoll/^ol + + laall/^al) " ^laaPlAP + ^(KWPol + laJIAD' 

-^(l«ollA|-|«JI/3ol)'-^(l«.l' + 1/3.1')-^- 
By using the Cauchy-Schwartz inequality, we have 

< 6|a,||a,| - ^la^lW + ^[(l«oll/5ol + " (l«ollAI " MPolf] " ^(kJ' + lAl') " ^- 

The parametrization P-3P yields 

= 6 cos ^ cos 0' cos^ e cos^ 0' + ^ [sin^ 6 sin^ 0' cos 2{(p - - cos^ 9 - cos^ - ^. 

Since the Cauchy-Schwartz inequality becomes an equality when the two vectors are parallel, we must 
have 9 = 9' and ip = ip' . Then 

< 3cos2 6'-cos^6l- -. 

4 

The maximum value of the latter expression is which come from cos^ 9 = ^. 

Appendix III : 

Proving the (l3AE\i : 
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To prove p.lSp . we try to find all pure states that are orthogonal to all pure product states 0-4^ . To 
this aim, let us consider a pure state as 

2 

i,j=0 

and try to find the coefficients a-, such that {tp) be orthogonal to all pure product states ()I-4p . We 
note that 

(7IV') = floo (3 cos^ ip) + a,)^ (Se"**^! sin 99 cos ip) + a^^ (\/3e~*''2 cos (p) + a^^ (Se"**^! sin ip cos ip) 

+a^, (Se-^^-'i sin2 p) + a^^ {VSe''^^^ ) sin v?) + a,,,, {VSe''^^ cos v?) (1-5) 
+03, ( V3e-*(''i +^2 ) sin p) + e-^^-'i +^2 ) . 

For p = 0, we get 

3a„„ + \/3e-*'^2 (a„^ + aj + a^.e-'^'^^'^) = 
and arbitrariness of phases 6^ and 62 follows that 

99 = ^ yields 

3a,,e-^'^i + \/3e-^(^i+'^2)(a^^ + a^i) = 
and, as before, arbitrariness of phases 5-^ and 62 follows that 



Substituting the above results in (|I-5p . we have = Finally, it is concluded that the mutually 

orthogonal pure states 

K) = i=(|01)-|10)) , |V;,) = i=(|02)-|20)) , |V;3) = i=(|12)-|21)) (1-6) 
are orthogonal to all pure product states p-4p and the positive operator V in (j3.15|) is 

r = ^\^P,){i^A + li\A){i^2\ + \i^^){i^,\)■ 



Now, we consider another pure state as 



i,j=0 



Two-qutrit EWs 



25 



and try to find the coefficients b.^ such that |(^) be orthogonal to all pure product states \a*) \P) in 
P-4p . By similar arguments as above, we conclude that 

l¥'> = ^ (|00) + |11)-3|22)) 

and 
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